Measuring the Kondo effect in the Aharonov-Bohm interferometer 
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The conductance Q of an Aharonov-Bohm interferometer (ABI), with a strongly correlated quan- 
tum dot on one arm, is expressed in terms of the dot Green function, Gdd, the magnetic flux (p and 
the non-interacting parameters of the ABI. We show that one can extract Gdd from the observed 
oscillations of Q with (j>, for both closed and open ABI's. In the latter case, the phase shift /3 deduced 
from Q ~ A + B cos{(j} + /3) depends strongly on the ABI's parameters, and usually (3 ^ Tr/2. These 
parameters may also reduce the Kondo temperature, eliminating the Kondo behavior. 

PACS numbers: 73.21.-b, 73.23.-b, 71.27.-ha, 72.10.Fk 
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The recent observation of the Kondo effect in quan- 
tum dots (QD's), whose parameters can be tuned con- 
tinuously [1] , has been followed by much theoretical and 
experimental activity. For temperatures T below the 
Kondo temperature Tk, the spin of an electron local- 
ized on the QD is dynamically screened by the electrons 
in the Fermi sea, yielding a large conductance G through 
the QD, close to the unitary value 2e^/h, and a trans- 
mission phase a equal to 7r/2 [2,3]. A good tool to 
test these predictions involves embedding the strongly- 
correlated QD on one arm of an Aharonov-Bohm inter- 
ferometer (ABI). Indeed, such experiments were carried 
out for both a closed (two-terminal) ABI [4] and an open 
(multi-terminal) ABI [5]. Both experiments exhibited 
the Aharonov-Bohm oscillations with the normalized flux 
(j) = e^/hc. The former experiments exhibited the ex- 
pected "phase rigidity" , with Q an even function of 4> 
[6]. However, there has been no quantitative analysis of 
these data. The latter experiments attempted to measure 
the transmission phase, and found a variety of behaviors 
which were inconsistent with the expected value of 7r/2. 
As a result, Ji et al. [5] stated that "the full explanation 
of the Kondo effect may go beyond the framework of the 
Anderson model". Theoretical attempts to discuss re- 
lated issues have concentrated on the dot alone (when it 
is detached from the ABI) [3,7], or apphed various tech- 
niques [8-10] to the QD on simple models of the closed 
ABI. However, it has not been very clear how to make 
quantitative comparisons of theory and experiment. 

Most of the theoretical discussions of QD's concentrate 
on the retarded Green function for electrons with energy 
Lo on the QD, Gdd{^) (we ignore the spin index, since 
we assume no magnetic asymmetry). For a simple QD, 
connected to a broad electronic band, the T — Q trans- 
mission amplitude for electrons going through the QD is 
proportional to Gddi^p), where e^? is the Fermi energy 
(taken as zero below) [2,11]. However, measuring the 
conductance only yields \Gdd{^)\'^, with no information 
on the phase. The ABI experiments were thus intended 
to measure both the magnitude and the phase of Gdd, 
and compare with theory. In this paper we concentrate 



on the following question: given experimental data on the 
flux dependent conductance of the ABI, Q{4>): how can 
we deduce the "intrinsic" Green function G^d? An ear- 
lier paper [12] answered this question for non-interacting 
electrons on a simple model for a closed ABI, and made 
some speculations on the interacting case. Another paper 
[13] showed that for non-interacting electrons, the phase 
shift measured in the open ABI depends on details of the 
opening. Here we show that for strongly correlated elec- 
trons, Q{<j)) is much more sensitive to the specific details 
of the ABI's. We give explicit instructions for extracting 
Gdd from the measured 5(0), and show that the opening 
has much stronger effects in the Kondo regime, possibly 
explaining the puzzling experiments [5]. 



Our qualitative results should apply for a large class 
of ABI's. For simplicity, we demonstrate them for the 
specific Anderson model shown in Fig. 1, which cap- 
tures the important ingredients. The conductance Q is 
measured between the two leads which are attached to 
sites "L" and "R" on the ABI ring. The QD (denoted 
"D") is connected to L (R) via ni (n^) sites. The lower 
"reference" branch contains sites. Except for the QD, 
we use a tight binding model, with the real hopping ma- 
trix elements as indicated in the figure. Site energies 
are e;, and eo on the respective branches, e^, eR on 
sites L and R and zero on the leads. Using gauge in- 
variance, we introduce the normalized flux (j) as a phase 
factor in Joi = Jlo ~ j/e"*. The Hamiltonian on the 
dot is Tid = Erf X^cr "■rfc + Uud^ridi, with obvious nota- 
tions. Here we assume that the transport is dominated 
by the level on the QD. We also assume that U is 
very large, and ignore the resonance at 2ed + U . Figure 
1 generalizes earlier models [8,9], by adding the internal 
structure on the links between D, L and R. Some such 
structure always exists in experiments, and may have im- 
portant effects on the observed conductance (see below) . 
For the open ABI, each dashed line represents an addi- 
tional lead, with a hopping matrix element —Jx on its 
first bond [13]. 
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FIG. 1. The model for the ABI, with m = rir = 3, no = 5. 

For T <C Tk, it is sufficient to calculate ^ at T = 0. 
Irrespective of the above details, one has 

g = ^4rL{o)TR{o)\GLRiio = o)\\ (1) 

where Glr{co) is the Fourier transform of the retarded 
Green function involving a creation of an electron at site 
R and an annihilation of an electron at site L at time t 
later [14]. Also, Tl,r{<^) = -3S2,^^(w), where E?,^^ is 
the self-energy generated at sites L or R due to the leads. 
Since interactions exist only on the dot, one can use the 
equations of motion to express all the retarded Green 
functions (j„/3 in terms of Gdd- If .9q/3 denotes the Green 
functions of the whole network without the QD and the 
two bonds connected to it, then it is straightforward to 
obtain [15] the relation 

GlR = QLR + ^ gLaJaDJD(}90RGdd, (2) 
a/3 

where in our model the only non-zero Ja^s are Jdi = 
JiD iie"^ and Jdr = Jrd = jr (sec Fig. 1). 
For the non-interacting case, one has 

Gdd = l/[w-ed-EoH], (3) 

with the non-interacting self-energy 

So(w) = ^ JDaga0{(^)J0D = Sed{u)) - iAo{u)). (4) 

a/3 

Since the phase (j> is only contained in Jui = Ji^, 
the complex matrix ga/s is independent of (j) and obeys 
5a/3 =g0a- Therefore, 

^0 (t^) = if 511 + JrdNN + '^3l3r9lN COS (j) (5) 

is an even function of (j). Except for the very special case 

ni = rir = riQ = (discussed e.g. in Refs. [8,9]), giN is 
not a real number, and therefore both Sed{w) and Ao(a;) 
oscillate with (f), 

Sed = ai + bi cos (j); Aq = 02 + 62 cos (/>. (6) 

Using these expressions, Eq. (2) can be written as 

GlR = gLRGdd([Gdd]~^ + So -h ^ JDaWa^Jpo^-, (7) 

a/3 



with Waf} = gaRgLpI gLR. — gaf3- In Calculating the nec- 
essary ga/3^s, it is convenient to first eliminate all the 
leads, replacing them by site self-energies. For example, 
at site L we replace the bare Green function l/{uj — e^) 
by — cl — S^(w)]. For our one-dimensional leads 
one has = -e'l^l" j£/J, where w = -2Jcos(q'a) is 

the energy of the electron in the band of the leads (we 
assume all the leads to have the same J and the same 
lattice constant a) [15]. The results remain valid also 
for more complex leads, as long as one uses the appro- 
priate band-generated self-energies. For the remaining 
A?' = -|- -|- no -|- 2 sites of the ring (without the QD), 
gaf} is then the inverse of a tri-diagonal N x N matrix, 
■Mai3 = Ja/3 + [<^ ~ Ca " ^a('-^)]^a/3- For such a matrix, 
it turns out that one always has win = 0. Thus, 

g = GreflGddflGd^ + i:o + x + ye-''^\\ (8) 

where x = jf 'Wn + jf-WNN-, y = jijrWNi and the "refer- 
ence" conductance g„/ = (2e2//i)4rz,(0)rij(0)l5Lfi|^ de- 
pend only on the parameters of the non-interacting parts 
of the ABI, and not on the QD parameters and U. 

At T = the electrons must obey the Fermi liquid re- 
lations [16]. Specifically, at the Fermi energy one expects 

S[Grfi](a; = 0) = Ao(a; = 0). (9) 

This condition should hold for any network in which 

the QD is embedded, and is therefore true for both the 
closed and the open ABI. In the limit of a very large 
negative e^, when (rid) — >^ 1) one also expects that 
^[G^^]{u; = 0) ^ 0. For the simple QD with two leads, 
one has So(a;) = S2,(w) -|- T,'^{lu), and in the symmet- 
ric case Tl = Tr = Ao/2 the conductance reaches its 
unitary limit 2e'^ /h [Eq. (1), with Glr Gdd]- In this 
limit the phase of Gdd becomes n/2 [2]. 

We now discuss the closed ABI. In this case, wn , wjvat 
and wni and therefore also x and y are all real numbers. 
Using the Fermi liquid result (9) and Eq. (6), Eq. (8) 
becomes 

gciosed _ , X _ (C + + n COS (t)f + rl sin^ 

Qref --^^^P'- C2 + (l+r-rfCOS</.)2 ' ^'"^ 

with the dimensionless function (^{(p) — 5R[G'^J]/a2 and 
constants Tq = (ai -I- x)/a2, rb — (6i -I- y)/a2, ry = y/a2 
and rd = &2/a2. For the non-interacting case, C = 
[—€d — Sed{0)]/a2, and Eq. (10) generalizes the results 
of Ref. [12]. An example is shown on the LHS of Fig. 2. 
In the strongly correlated case and in the unitary limit, 
C = and Qclosed ^ Qre}^o{4')- All the features in the 
(^-dependence of J^o arise only due to the non-interacting 
parts of the ABI. Usually, Eq. (10) contains many har- 
monics. Except in special cases [8], it is not dominated 
by the second harmonic, and the period of To{(t>) is not 
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simply doubled. An example of this dependence is seen 
(for large negative e^) on the RHS of Fig. 2: except 
for the minima at = and tt, the maxima are not at 
7r/2. Experimentally, one knows that one has reached 
this limit once the function Qciosed{<P) no longer changes 
with the gate voltage which governs e^. The reference 
conductance Qref can be measured by disconnecting the 
QD, i.e. setting ji — jr = 0. Alternatively, Gref can 
be absorbed in the scales of the parameters in the nu- 
merator of Eq. (10). Having determined Grefj one can 
determine the four real parameters ra, ry and by 
a fit to J-o{(j)) (In practice, one only needs four values of 
the function) [17]. Having found these parameters, one 
can now move away from the unitary limit, and measure 
Gciosed = Gref^{(t>)- The uukuowu fuuctiou ({(p) can now 
be found from the quadratic equation 
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rg + n cos ( 



(1 + Td COS( 



0. (11) 



The solution should be chosen so that it decreases to zero 
at large negative and increases linearly with large pos- 
itive ed- Having found the solution, the phase a of Gdd 
is then defined via 



cot a — 



jR[G-J](^ = 0) 
Ao(c^ = 0) 



C 



I + rd cos ( 



(12) 



This phase, or equivalently 5ft[G'^J], are the quantities 
obtained from theories. 





FIG. 2. Conductance (in units of 2e^ /h) through the closed 
ABI versus the normalized flux (f} and the energy of the state 
on the dot ed (the gate voltage), without (LHS) and with in- 
teractions (RHS). m — rir = 2, no = 3, Jl = Jr — Jd — 1, 
ji = jr = -2, ii = ir = .4, e; = = eo = --3, el = eji = 0. 
All energies are in units of J. 

For demonstrating the qualitative dependence of G on 
(j) and on the other parameters, we have used an ap- 
proximate analytic solution of the equations of motion, 
truncated via decoupling of higher order Green functions 
[15]. In the limits T = lj = Q and C/ — > oo, this solution 
assumes the simple analytic form 



cot a 



5n + z(^z + \J z"^ + '^'^(i ~ f^"-)) 



(13) 



where z represents the value at = of the non- 
interacting ratio 

z(c^) = [oj-ed~ 5ed{u)]/[2^a{uo)], (14) 

while 5n is related to the electron occupation on the dot 
via {ud) = 2(1 — 5n) (which should be determined self- 
consistently) . In practice, 5n varies smoothly between 
1/2 (at z 1) and 1 (at z ^ —1), and the results of 
calculations are not very sensitive to the details of this 
variation. Equation (13) interpolates between /3 = 7r/2 at 
z ^ oo and f3 = n for z — oo. In the latter limit, Gdd 
approaches the non-interacting form (3). Using Eq. (13) 
in Eq. (10) for a specific set of parameters yields the RHS 
of Fig. 2. One clearly sees the transition from the non- 
interacting behavior at large positive ed (compare with 
the LHS) to the unitary limit at large negative ed- For 
different sets of parameters one reproduces qualitatively 
all the earlier results, including the Fano-Kondo effect 
[8]. We have used these results to imitate real experi- 
mental "data" , and were able to use the above algorithm 
to extract cot a as in Eq. (12). 

Note that the above analysis yields Gdd for the QD on 
the ABI, where this function (and thus also the phase a) 
depends explicitly on the ffux (/), via z. At T = w = 0, 
we expect a to depend only on the ratio z also for other 
theories. In our case, z can be extracted from the exper- 
imental data via 

z = -[id + ra + {n - Ty) coS(/)]/(l + Tdcoscj)), (15) 

where Id — {ed — x) j a-i is just a shifted rescaled gate volt- 
age. Having deduced the dependence of both z and a on 
(/), a parametric plot can yield a versus z, for compari- 
son with single dot calculations. Alternatively, one can 
experimentally study the results as function of the cou- 
pling to the reference branch, ii and ir- Extrapolation to 
ij, V ^ would give the dependence of a on So(0) for 
the upper branch alone. However, So('-^) still depends on 
the finite chains connecting D with L and R [18]. 

We now turn to the open ABI, with Jx / 0. Equation 
(8) remains correct, but now x and y become complex. 
Interestingly, Eq. (5) still holds, and Eq is still an even 
function of (j). In the unitary limit, Gi4') hS'S the exact 
form 

A + B cos(,/> + P) + C cos(2,/> + 7) 

yopen-^yref (1 + COS (/.)2 ' 

and we need six parameters to fit it. Note that all the 
ABI parameters (including Gref) now also depend on Jx- 
The two lower curves in the left panel of Fig. 3 show 
results in this limit. Note that the graphs are not sinu- 
soidal, mainly due to the second term in the numerator 
and to the denominator in Eq. (16). Since one remains 
close to the Kondo resonance, the denominator contin- 
ues to be important, modifying the 2-slit-like numerator. 
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The asymmetric shape of each oscillation seems similar 
to that reported in Ref. [5]. The other curves in the 
same panel were derived using Eq. (13). Again, one 
observes the crossover to the non-interacting sinusoidal 
shape at large positive e^. To extract a "transmission 
phase" from these curves, one can e.g. follow the max- 
ima as function of Cd, or enforce a fit to the two-slit for- 
mula Qopen ~ ^ + Bcos{(j) + j3). Since now there is no 
well-defined zero to (j), one can only deduce the relative 
change in the phase (3. Setting /3 = at ^ — oo, the 
RHS of Fig. 3 shows this relative phase versus e^. For 
the parameters we used, the total change is about O.Stt, 
far away from the expected change in a, equal to 7r/2. 
The actual values depend on details of the ABI. This may 
explain the non-trivial values of the phases observed in 
Ref. [5]: they result from the experimental setup, and 
not from a breakdown of the Anderson theory. 
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FIG. 3. LHS: Conductance through the open ABI versus 
</>, at ed = ( — 1.5, —1, ^.5, 0, .5, 1, 1.5) J, with interactions. 
Graphs are shifted up with increasing ed- Parameters are the 
same as in Fig. 2, but with Jx = .5. RHS: The "measured" 
phase shift /3 (in units of tt) versus ea- 

Finally, a few words about non-zero T or w. Generally, 
T and a; enter into Gm similarly. In the approximate 
solution of Ref. [15], one ends up with a competition be- 
tween the variable z{u}) of Eq. (14) and log(£'/T) or 
\og{D/uj), where 2D = AJ is the width of the band in 
the leads. This competition yields estimates of Tk, 

\og{TK/D) = TT[ea + .5ed(0)]/Ao(0). (17) 

Although more accurate theories end up with different 
expressions, all of them end up with a strong dependence 
on the ratio which appears on the RHS. In our case, this 
ratio oscillates strongly with 0, opening the possibility 
that for different fluxes the QD is below or above Tk- 
We emphasize the appearance of 5ed in the numerator, 
ignored in some papers. 

At non-zero T, the "intrinsic" phase of the QD is ex- 
pected to start at for large negative [where T > 
TRi^d)], then grow to 7r/2 for intermediate negative e^'s 
(the unitary region), and finally grow to tt at positive 
Crf [3]. As mentioned, both Sea and Aq depend on the 
opening parameter Jx ■ Using the approximation of Ref. 
[15] also for T > 0, we found that large values of Jx may 
completely eliminate the intermediate plateau in a, and 
give a direct increase of a from to tt. Unlike the non- 



interacting case [13], where changing Jx only slightly 
modified the quantitative shape of the function j3{ed), 
the effects here are qualitative: opening may lower Tk 
and completely eliminate the observability of the Kondo 
behavior. Again, this could have happened in Ref. [5]. 
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